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We study the influence of quenched disorder on quantum phase transitions in systems with over- 
damped dynamics. For Ising order parameter symmetry disorder destroys the sharp phase transition 
by rounding because a static order parameter can develop on rare spatial regions. This leads to an 
exponential dependence of the order parameter on the coupling constant. At finite temperatures 
the static order on the rare regions is destroyed. This restores the phase transition and leads to 
a double-exponential relation between critical temperature and coupling strength. We discuss the 
behavior based on Lifshitz-tail arguments and illustrate the results by simulations of a model system. 



The influence of quenched disorder on phase tran- 
sitions is an important problem in condensed matter 
physics. Initially it was suspected that disorder destroys 
any critical point However, it was soon found that 
classical continuous phase transitions generically remain 
sharp in the presence of weak disorder. If a clean critical 
fixed point (FP) fulfills the Harris criterion |^ v > 2/d, 
where u is the correlation length exponent and d is the 
spatial dimensionality, weak disorder is renormalization 
group irrelevant. The system becomes asymptotically ho- 
mogeneous at large length scales. Even if the Harris cri- 
terion is violated, the transition will generically be sharp. 
In this case the inhomogeneities either remain finite at 
all length scales, leading to a finite disorder FP, or the 
relative magnitude of the disorder diverges under coarse 
graining corresponding to an infinite-randomness FP. A 
prominent example of the latter occurs in the McCoy- 
Wu model a disordered 2D Ising-model in which the 
disorder is perfectly correlated in one dimension. These 
correlations increase the effects of the disorder. 

An important aspect of phase transitions in disordered 
systems are the Griffiths phenomena Q . They are caused 
by large spatial regions that are devoid of any impurities 
and can be locally in the ordered phase even if the bulk 
system is in the disordered phase. Since these regions are 
of finite size, no true static order develops. The fluctua- 
tions of these regions are very slow because they require 
changing the order parameter in a large volume. Griffiths 
[Ql showed that this leads to a singular free energy. In 
generic classical systems this is a weak effect, since the 
singularity is only an essential one. An exception is the 
McCoy- Wu model Here, the disorder correlations 
lead to stronger effects, with the average susceptibility 
diverging in a finite region of the phase diagram. 

Recently disorder effects have gained a lot of atten- 
tion in the context of quantum phase transitions |^ . At 
these transitions order-parameter fluctuations in space 
and time have to be considered. Quenched disorder is 
time-independent, it is thus always correlated in one of 
the relevant dimensions making disorder effects at quan- 
tum phase transitions generically stronger than at clas- 
sical transitions. A prototypical model is the random 
quantum Ising fcrromagnet. Its quantum phase tran- 



sition in ID §, |, I and 2D §, |T§ is controlled by 
an infinite-randomness FP with activated rather than 
power-law dynamical scaling. The Griffiths singularities 
are enhanced, too: Several observables including the av- 
erage susceptibility display power-law singularities with 
continuously varying exponents over a finite region of 
the disordered phase. Similar phenomena have also been 
found in quantum Ising spin glasses Jll| , |l^ , [Tst . 

The systems in which infinite-randomness FPs and 
quantum Griffiths phenomena have been shown unam- 
biguously all have undamped dynamics (a dynamical ex- 
ponent z = 1 in the corresponding clean system). How- 
ever, in itinerant electronic systems the order parame- 
ter couples to fermionic modes leading to overdamped 
dynamics (clean dynamical exponent z > 1). A proto- 
type example is the itinerant quantum antifcrromagnetic 
phase transition. The conventional pcrturbative renor- 



malization group for this transition |14| yields a finite- 
disorder FP. However, by taking into account the effects 
of rare regions it was shown |15|| that this FP is unsta- 
ble, and the renormalization group flow is towards large 
disorder. The meaning of this runaway flow is presently 
not fully understood. In addition to the transition itself, 
quantum Griffiths phenomena in itinerant systems have 
also attracted much attention since they are of potential 
importance for a variety of heavy- fermion systems. It 
has been suggested that overdamped systems show 
quantum Griffiths phenomena very similar to that of un- 
damped systems. However, recently it has been argued 
|l7t that for Ising symmetry the overdamping prevents 
the rare regions from tunneling leading to superpara- 
magnetic rather than quantum Griffiths behavior. 

In this Letter we reconsider the important question 
of disorder and rare regions at quantum phase transi- 
tions. We show that for Ising order parameter symmetry 
and Landau overdamped dynamics, the sharp quantum 
phase transition is destroyed by rounding, because static 
order can develop on isolated rare regions. Our results 
can be summarized as follows: The relation between the 
order parameter m and the distance t from the clean 
critical point is exponential for small m. The precise 
form depends on the disorder distribution. For a Gaus- 
sian coupling constant distribution, m is finite for all t. 
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Asymptotically it behaves as 

log(TO) - -t^-'i/^ for i ^ cx) . (la) 

Here 4> is the finite-size scaling shift exponent of the clean 
system. If the disorder is of dilution type a finite order 
parameter starts to develop at the clean transition, i.e., 
for t < For small \t\ it behaves as 

log(m) \t\-'^/<l' fort^O-. (lb) 

At finite temperatures the static order on the rare re- 
gions is destroyed, and a finite interaction between them 
is required for long-range order. This restores a sharp 
phase transition. The dependence of the critical temper- 
ature on t is double-exponential. For small Tc we obtain 

log(-alogTj - f~'^''>' Gaussian, (2a) 
log(-alogr,) ~ \t\-'^''l' dilution, (2b) 

Moreover, in finite-size samples strong sample-to-sample 
fluctuations occur if the number of rare regions displaying 
static order becomes of order one. Asymptotically for 
large linear system size L, the coupling strength at 
which these fluctuations start behaves like 

tL - (logi)i/(2-d/0) Gaussian, (3a) 
\tL\ - (logi)-"*/'' dilution. (3b) 

Thus, finite size efi^ects are suppressed only logarith- 
mically. In the remainder of this Letter, we sketch the 
derivation of these results, illustrate them by numerical 
results from a model system and discuss their relevance 
for experiments as well as simulations. 

For definiteness we consider the antiferromagnetic 
quantum phase transition of itinerant electrons with Ising 
spin symmetry. The Landau-Ginzburg- Wilson free en- 
ergy functional of the clean transition reads ||lj, 

S = J dx dy m{x)T{x,y) m{y) + u J dx m^{x) (4) 

Here x = (x, r) comprises position x and imaginary 

time T, and J dx = J dx J^^^ dr. T{x,y) is the bare 
two-point vertex, whose Fourier transform is r(q, a;„) — 
{t + + \LJn\) ■ Here t — {g — gc)/gc is the distance of the 
coupling constant g from the critical point. The dynam- 
ical part of F is proportional to |w„| reflecting the over- 
damping of the dynamics (undamped dynamics leads to 
w^). Quenched disorder is introduced by making t a ran- 
dom function of position, t — *■ t + (5i(x). We consider two 
different disorder distributions. The flrst is a Gaussian 
distribution with zero mean and a correlation function 
((5t(x)(5<(y)) = A^5(x — y). The second disorder distri- 
bution is of dilution type, (5<(x) = everywhere except 
on randomly distributed finite-size islands (impurities) of 
spatial density c where (5i(x) = W > Q. 



In the presence of disorder there are rare large spatial 
regions which are locally in the ordered phase while the 
bulk system is not. For the Gaussian disorder distribu- 
tion which is unbounded, these regions exist for all t. For 
the dilution case they appear below the transition of the 
clean system, i.e. for t < 0. At zero temperature a sin- 
gle such region is equivalent to a classical Ising model in 
a rod-like geometry. It is finite in the d space dimen- 
sions but infinite in imaginary time. If the interaction in 
the time direction is short-ranged, as is the case for un- 
damped dynamics, true static order cannot develop on 
such a rare region. Instead, the order parameter displays 
slow fluctuations leading to quantum Griffiths phenom- 
ena. This is drastically different in a system with over- 
damped dynamics. The linear frequency dependence in 
F is equivalent to a long-range interaction in time of the 
form (r — t')~^. ID Ising models with 1/r^ interaction 
are known to develop long-range order at finite temper- 
atures |po| . Thus, in our quantum system, true static 
order develops on those rare regions which are locally 
in the ordered phase. In agreement with Ref. |l^ we 
therefore do not find the usual quantum Griffiths phe- 
nomena. Once static order has developed on a few iso- 
lated rare regions, an infinitesimally small interaction or 
an infinitesimally small symmetry-breaking field are suf- 
ficient to align them. Consequently, a macroscopic order 
parameter arises. 

We now use Lifshitz tail arguments |2l] to derive the 
leading thermodynamic behavior for small order param- 
eter m. In the dilution case, the probability w to find a 
region of linear size devoid of any impurities {5t = 0) 
is given by w ~ exp(— cL^) (up to pre-exponential fac- 
tors). Such a rare region develops static order at a some 
tciLn) < 0. Finite size scaling yields \tc{Lj^)\ ~ 
where (j) is the finite-size scaling shift exponent of the 
clean system Thus, the probability for finding a 

rare region which becomes critical at tc is given by 

w{tc) ^ expi^B \tc\^'^/^) fort^O-. (5a) 

For a Gaussian distribution similar arguments give 

w{tc) ~ exp(-B t^-''/'^) for t ^ oo . (5b) 

Here B and B are constants. The total order parameter 
m is obtained by integrating over all rare regions which 
are ordered at t, i.e., all rare regions having tc > t. This 
leads to eqs. (|l|). Note that the functional dependence on 
t of the order parameter on a given island is of power-law 
type and thus only influences the pre-exponential factors. 

At finite temperatures the static order on the rare re- 
gions is destroyed, and a finite interaction of the order of 
the temperature is necessary to align them. This means 
a sharp phase transition is recovered. To estimate the 
transition temperature we note that the interaction be- 
tween two rare regions depends exponentially on their 
spatial distance r, Ei„t ~ exp(— r/^o)i where is the 
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bulk correlation length. The typical distance r itself de- 



pends exponentially on t via r 



w 



-l/d 



The leading 



dependence of the critical temperature on t is thus 

^ cxp(-r/eo) ^ cxp(-a nj-'^^/^o) , (6) 

where a is a constant. Inserting eqs. for w leads to eqs. 
(^. Above this exponentially small critical temperature 
the rare regions essentially behave classically. 

We now turn to finite-size effects at zero temperature. 
Since the total order parameter is the sum of contribu- 
tions of many independent islands, finite size-effects in 
a macroscopic sample are governed by the central limit 
theorem. However, for t — > oo (Gaussian distribution) or 
t — > 0— (dilution) very large and thus very rare regions 
are responsible for the order parameter. The number N 
of rare regions which start to order at tc in a sample of 
size L behaves like TV ~ L'^ w{tc)- When N becomes 
of order one, strong sample-to-sample fluctuations arise. 
Using eqs. (||) for w{tc), this leads to eqs. (^). 

To illustrate the rounding of the phase transition we 
now present numerical data for a model system, viz., a 
disordered classical Ising model with two space dimen- 
sions and one time-like dimension. The disorder is of di- 
lution type and completely correlated in time direction. 
The interaction is short-ranged in the space directions 
but infinite-ranged in the time-like direction. This sim- 
plification retains the crucial property of static order on 
the rare regions but allows us to treat system sizes large 
enough for the investigation of exponentially rare events. 
The Hamiltonian of the model reads 



H = - 



Lr 



Here x and r are the space and time coordinates, respec- 
tively. Lt is the system size in time direction and (x, y) 
denotes pairs of nearest neighbors. Jx is a quenched 
binary random variable with the distribution P(J) = 
(1 — c) S{J —l) + c S{J). In this classical model Lr takes 
the role of the inverse temperature in the corresponding 
quantum system and the classical temperature takes the 
role of the coupling constant g. Because the interaction 
is infinite-ranged in time, the time-like dimension can be 
treated in mean- field theory. For Lr ^ oo, this leads to 
a set of coupled mean-field equations (one for each x) 



mx = tanh f3 [J^fn^ 



y(x) 



h] 



(8) 



where we have added a very small symmetry-breaking 
magnetic field, h = 10^^. These equations are solved 
numerically in a self-consistency cycle. 

Fig. |l| shows the total magnetization and the suscepti- 
bility (which corresponds to the inverse energy gap of the 
quantum system) as functions of temperature for linear 
size L — 100 and impurity concentration c = 0.2. The 
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FIG. 1: Magnetization and susceptibility {L — 100, c = 0.2). 
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FIG. 2: Magnetization and susceptibility close to the seeming 
transition for different system sizes. 



data are averages over 200 disorder realizations [g4j. At 
a first glance these data suggest a sharp phase transition 
close to r = 4.88. However, a closer investigation, Fig. ||, 
shows that the singularities are rounded. If this round- 
ing was a conventional finite-size effect the magnetization 
curve should become sharper with increasing L and the 
susceptibility peak should diverge. This is not the case 
here, instead the transition remains rounded for L — > oo. 

For comparison with the analytical results. Fig. ^ 
shows the logarithm of the average magnetization as a 
function of 1/(T° - T) where T° = 5 is the critical tem- 
perature of the clean system (c = 0). The data follow 
eq. ( p^ ) over almost four orders of magnitude in m with 
the expected shift exponent oi (p — 2. This figure also 
shows "typical", i.e., logarithmically averaged magneti- 
zation data for different system sizes. Deviations between 
the typical and the average values (which are essentially 
size-independent) reflect strong sample-to-sample fluctu- 
ations. The data show that the onset Tl of these fluctua- 
tions shifts to larger temperatures with increasing system 
size. A more detailed analysis shows that = (T^ — T°) 
follows eq. ( |3b| ) in good approximation. 

In the remaining paragraphs we discuss the general- 
ity of the results and their consequences for experiments 
and simulations. The rounding of the transition is due 
to the fact that at zero temperature static order can de- 
velop on a finite spatial region in a Landau-damped sys- 
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FIG. 3: Log(m) as a function of distance from tiie clean crit- 
ical point. The solid line is a fit of the average magnetization 
to eq. ( ]lb[ ) with (f> = 2. The logarithmically averaged data 
show the onset of the sample-to-sample fluctuations. 



tern. We thus expect all quantum phase transitions with 
discrete order parameter symmetry and overdamped dy- 
namics (i.e., a fow frequency dependence ~ |a;| or slower 
in the Gaussian propagator) to be rounded by disorder. 
Systems with continuous symmetry behave differently, ft 
is known that classical ID XY and Heisenberg sys- 
tems develop long-range order at finite T only if the in- 
teraction falls off more slowly than I/r^. Consequently, 
in a quantum system at zero temperature static order on 
a rare region only develops if the frequency dependence 
in the Gaussian propagator is slower than |u;|. The itin- 
erant ferromagnetic or antiferromagnetic quantum phase 
transitions with XY or Heisenberg symmetries will thus 
remain sharp in the presence of disorder . 

The rounding of the quantum phase transition leads to 
an unusual phenomenology in experiments. Data taken 
at larger order parameter and not too low temperature 
do not resolve the exponentially small order parameter 
tail but probe the rounded transition. However, with 
increasing precision and decreasing temperature the ap- 
parent critical point shifts further and further towards 
the disordered phase, accompanied by the onset of strong 
sample-sample fluctuations. Similar effects will also oc- 
cur in simulations, fn a recent Monte-Carlo simulation 
of a related model ||2^ a sharp transition was found con- 
trolled by an infinite randomness fixed point. However, in 
this simulation the linear system sizes are below L = 30, 
therefore it likely probes the rounded transition rather 
than the exponentially small magnetization tail. 

In this Letter we have concentrated on the behavior of 
the system for very small order parameter values, i.e., in 
the "tail region" of the rounded transition. The above ar- 
guments suggest that the properties of the rounded tran- 
sition itself are experimentally also very important be- 
cause they control the behavior in a wide pre-asymptotic 
region. However, these properties are likely to be non- 
universal, so far they are not very well understood. 
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